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I review work on diffractive vector meson production in photon— proton collisions at high 
energy and large momentum transfer, accompanied by proton dissociation and a large 
rapidity gap. This process provides a test of the high energy scattering dynamics, but is 
also sensitive to the details of the treatment of the vector meson vertex. 

The emphasis is on the description of the process by a solution of the non-forward 
BFKL equation, i.e. the equation describing the evolution of scattering amplitudes in the 
high-energy limit of QCD. The formation of the vector meson and the non-perturbative 
modeling needed is also briefly discussed. 
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1. Introduction 

High energy reactions producing mesons, and in particular diflractive processes, 
are advantageous for studying various aspects of high energy scattering processes, 
partly because such processes are experimentally clean, and partly because many 
processes offer unique types of reactions. 

This brief review is about a certain class of meson production, namely diffractive 
vector meson production in photon-proton collisions, where the color singlet t- 
channel exchange carries large momentum transfer, i.e., jp — > VX (see Fig. [Q. 
In such processes there is a large rapidity gap centrally in the detector and the 
proton dissociates. The vector meson is on one side of the gap and on the other 
side of the gap there will be a jet, balancing the transverse momentum. This jet 
is not measured in experiments, but the cross section is obtained simply from the 
kinematics of the vector meson — a very nice experimental signal. The cross section 
is rather small, however, and the available HERA datcPEl 

is therefore very recent. 

The mesons are identified experimentally by their decays into charged particles, 
e.g., p — ► 7r + 7r~ or J/tjj — > l + l~ , which provides a means of measuring the relative 
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Fig. 1. Diffractive photoproduction ■yp — > VX of a vector meson at large momentum transfer. 

polarization by examining the angular distributions of the decay products. a 

This process is quite different from common diffractive processes where the pro- 
ton survives intact and the rapidity gap is forward in rapidity. In general for pho- 
toproduction there can be two large momentum scales in the process — the squared 
momentum transfer t, and the vector meson mass My. In normal diffractive pro- 
cesses, t = 0, and only large-mass meson production can be treated perturbatively. 
But if —t 3> Aq CD perturbation theory can be used even for light mesons. 

The i-channel exchange must be in a color singlet state and is often referred to 
as hard color singlet exchange, or hard pomeron exchange. The lowest order contri- 
bution is therefore the exchange of two gluons. At high center of mass energy y/s, 
however, the amplitude is dominated by terms proportional to [a s ln(— s/t)] n and a 
fixed order calculation is not adequate. Rather, one needs to resum these logarithmic 
terms to all orders n. This is done by solving the Balitsky-Fadin-Kuraev-Lipatov 
(BFKL) equationSl, which resums the above logarithms in the leading logarithmic 
(LL) approximation. 

The BFKL equation describes the evolution of QCD amplitudes in the high 
energy limit s 3> — t and gives a new type of dynamics compared to fixed-order 
calculations. The resummation corresponds to considering diagrams with exchange 
of two gluons to all orders and keeping only those pieces containing a large logarithm. 
The main effects are that the two gluons are reggeized due to virtual corrections, 
which means that the propagators acquire a Regge-like factor with a power of — s/t, 
and that any number of gluons is exchanged between the two gluons such that a 
gluon ladder is formed. This gives a cross section which goes as a power of the 
energy, like in Regge theory. See e.g. Refs. |SJ 0for more on BFKL physics. 

The status of the BFKL resummation framework is not yet clear because of 
uncertainties related to experimental limitations, small predicted effects at present 
accelerators, NLL corrections, and so on. Many of the suggested signatures of BFKL 
resummation, such as the rise of F2 at small x, can also be described by the tra- 
ditional DGLAP equations. This makes it important to investigate possible signals 

a The background to p production is therefore the prompt production process -yp — > Tr+Ti^X, 
which could in principle be computed with the methods reported here by replacing the distribution 
amplitude of the p meson by the generalized 7r+7r — distribution amplitude . 
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for the behavior predicted by BFKL. In Ref. 7J vector meson production at large 
t was suggested to be one suitable process. Anot her large-t process that has been 
advocated is large rapidity gaps between jetsPEI _____ 

This review is organized around our calculation J 1 ^ 1 1 ^ 1 of vector meson produc- 
tion using BFKL resummation in the LL approximation. We will make use of the 
non- forwa rd v ersion of the equation, which has a complete analytical solution due 
to Lipatov.^ This solution, and our calculations, will be discussed in Section __ 

But there is of course also a lot of interest in the vector meson production 
process in itself. The experimental situation is not well understood from theory. 
The measured differential cross section da/dt for light meson production shows a 
independence which goes roughly as 1/t 3 , and measurements of spin density matrix 
elements indicate that these mesons are predominantly transversely polarized; only 
10-20% of the cross section comes from longitudinally pola rized mes ons. This is 
in contrast to the lowest order perturbative QCD calculation^2D____ that predict 
da/dt ~ 1/t 3 for longitudinal mesons and da/dt ~ 1/t 4 for transverse mesons — 
production of transversely polarized mesons is suppressed by an extra power of 
t\ I will come back to these problems below and discuss various approaches to 
their solution. In particular I will show that our BFKL calculations give a natural 
explanation of these hierarchies. 



2. Review of the theory for vector meson production 

The process jp — > VX t ake s place by exchange of the hard pomeron as shown in 
Fig. Q] It has been showrP^ that at large momentum transfer the pomeron couples 
predominantly to individual partons in the proton. One may therefore write the 
cross section as a product of the parton level cross section and the normal parton 
distribution functions of the proton, 



f 



where G(xj,t) and qf(xj,t) are the gluon and quark distribution functions respec- 
tively. The struck parton in the proton initiates a jet in the proton hemisphere, and 
carries the fraction Xj of the longitudinal momentum of the incoming proton. We 
are working in the high energy regime, were the coupling of the pomeron to a quark 
is identical to the coupling to a gluon, apart from the appropriate color factor. This 
is the reason for the factor in brackets on the right hand side in Eq. JIJ. So, in the 
following, all attention will be focused on the quasi-elastic subprocess — > Vq. 



2.1. Factorization 

The first computation^— ^ of vector meson photoproduction at arbitrary momentum 
transfer were done considering the Born graph for the pomeron, representing it 
by two gluons in a color singlet state. The results were given in the impact factor 
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(a) (b) 

Fig. 2. (a) AB — > CD scattering process in impact factor representation. The upper and lower 
blocks represent the impact factors for the corresponding transitions with two gluons attached, 
and the middle block marked HCS (hard color singlet) represents the Green function for the two 
gluons. For pomeron exchange, this is the BFKL Green function, (b) The 7 — > V impact factor. 



representation of Cheng and WuM^ This means that one writes the amplitude for 
the high energy process AB — > CD on the form 

M AB^CD [q) = J rf 2 fc fc2 J_ ^ C (fe, q) <P B ^ D (k, q), (2) 

where <P A ^ C and <P B ^ D are the impact factors for the upper and lower parts of 
the diagram, respectively. b That is, they are the impact factors for the processes 
A — > C and B — ► D with two gluons carrying transverse momenta k and q — k 
attached, the gluons being in an overall color singlet state. 

This factorization is schematically depicted in Fig. EH- It applies in the high- 
energy limit s 3> — t, where one can perform the integrals over the gluon plus and 
minus lightcone momenta, leaving only the integrals over the transverse components. 
The impact factorization simplifies calculations a great deal, because one can treat 
the process in three disjoint parts; the two impact factors up and down in the 
diagram, and the exchange mechanism in between. In the case of two-gluon exchange 
the middle block gives just delta functions connecting the gluon momenta and colors, 
and the denominator containing the propagators. In the BFKL case to be discussed 
later, the block contains the evolution of the two-gluon system into a gluon ladder, 
and will give a more complicated function. 



b For a derivation of this result and the expressions for the impact factors, see Ref. 1181 
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In Fig. 03, I show the lowest order diagram for the 7 — * V impact factor. The 
process is viewed in a frame where the photon has large lightcone momentum q + and 
fluctuates into a quark-antiquark pair before the hard scattering. This pair forms 
a color dipole which scatters by exchanging two gluons, and then forms a vector 
meson, which is possible because for pomeron exchange the vector meson has the 
same quantum numbers as the photon. The fluctuation of the photon into the dipole 
is long-lived, and the transverse size of the pair is frozen during the scattering. 

Also here there is a factorization formula which applies in the diffractive limit 
of negligible quark lightcone minus-momentum. The impact factor is given by 



& 7 (\)v(\'){k,q) ~ du d 2 r<P^^, ) (u,r)T(k,q;r,u)1 r l(x) (u 7 r), (3) 



where $V (S^) denotes the meson (photon) wave function (which will be dis- 
cussed below), A (A') denotes the polarization of the photon (meson), u is 
the quark lightcone plus-momentum fraction of the meson momentum, r is the 
transverse quark-antiquark separation, and T(k, q] r, u) is the hard amplitude 
of the QCD dipole scattering, which, up to a constant factor, equals J dl P° lc — 
e %qru A _ g-ifcr) '1 _ e —%(q-h)r\ rpj^g j m p ac ^ f ac tor for the q — > q transition is, in 

the high energy limit (where one uses the eikonal approximation), just a constant. 
The above results are given in a mixed representation in transverse momenta 
for the gluons and transverse vectors in coordinate space for the quarks. They can 
of course be written in the more familiar momentum space by a suitable Fourier 
transform, but the present form turns out to be quite useful for BFKL calculations. 

2.2. The photon and meson wave functions 

The photon and meson wave functions give the probability amplitudes for the split- 
tings 7 — ► qq and V — > qq. The photon wave function can be straightforwardly 
calculated in QED (see e.g. Refs.l51l2T|l and I do not reproduce the result here. d 

The meson wave function is more complicated. It contains non-perturbative 
information on the transition of the quark-antiquark pair into a vector meson; thus a 
model is needed. This model will usually involve the distribution amplitudes (DAs) of 
the meson. Roughly speaking, the DA 4>(u) describes the distribution of longitudinal 
momentum fraction u of the quark in the meson. In the early papers on light meson 
productio n 1 1 their form was not known, but one used either the a sym ptotic 
form (j>(u) — 6m (1 — u), or a form suggested by Chernyak and Zhitnitsky^ based 
on QCD sum rules. These calculations have invariably set the quark mass to zero. 
For heavy mesons, on the other hand, a non-relativistic quarkonium picture of the 
meson has been used, where the quark and the antiquark have the same longitudinal 

c There arc some subtleties regarding the pomeron coupling to colorful particles. This was addressed 
by Mueller and TanjCH See also the discussion in Refs. Ho1 and l2Ul 

d Though, when calculating the impact factor one might instead calculate the hard scattering 
amplitude for the whole ■ygg — > qq process directly and convolute it with the meson wave function. 
This is actually what is done in most of the literature cited here. 
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momentum, <j>(u) — 5(u— 1/2). In the BFKL calculation a ho be described below, 
we unify light and heavy mesons by keeping the quark masses throughout. We are 
therefore able to describe all mesons within the same framework. 

From quite general considerations one finds that the meson wave func- 
tion \Py above is related to the meson-to-vacuum transition matrix elements 
(V(p, A) | -0 ( — r /2)r tr ip(r /2)\Q) where -T* are Dirac-matrix structures from the Fierz 
decomposition of the amplitude. These matrix elements are precisely those that 
appear in the lightcone definitions of the D As and have been calculated from QCD 
sum rules with the inclusion of higher twistspSESI extending the results of Ref. [231 
I will not list the expressions here, suffice it to say that the distribution amplitudes 
of the vector mesons appear, and that they allow the calculation of the production 
amplitudes for all meson polarizations and quark helicities. 

2.3. Light meson calculations 

The hard scattering amplitude and impact factors were calc ulat ed quite some time 
ago in lowest order QCD by Ginzburg, Panfil, and SerbopH w ho obtained the 
amplitudes for production of longitudinally polarized light mesons in photon-photon 
and photon-quark collisions. Ginzbu rg a nd Ivanov^l added transversely polarized 
mesons to the calculation, and Ivanov*^ extended it to include BFKL effects. 

These results, however, have the features mentioned in the Introduction: the 
longitudinal mesons dominate the cross section for light mesons — because of the 
chiral even nature of the quark-photon coupling for massless quarks, the production 
of transversely polarized mesons vanishes as the quark mass tends to zero. 

A solution to this problem was suggested by Ivanov, Kirschner, Schafer, and 
Szymanowskij^l who argued that the photon splitting may also have a non- 
perturbative chiral odd component which would enhance the production of trans- 
verse mesons. This chiral odd photon wave function was modeled on a meson-like 
wave function, with a coupling constant arising from the quark condensate {qq). 
The coupling constant is proportional to (qq) and the magnetic susceptibility of the 
vacuum, which were taken from QCD sum rules estimates. The chiral symmetry 
breaking therefore occurs also when neglecting quark masses. 

This approach enhances the production of transverse mesons, such that it numer- 
ically dominates in the region of moderate momentum transfer that is experimen- 
tally observed, but longitudinal meson production still dominates for asymptotically 
large momentum transfers. However, the transverse meson cross section still scales 

°See Ref s . and 1 1 81 for more details. 

f The chirality of the qq state refers to the spinorial structure of the wave function — upon Fierz 
decomposing the amplitude, one gets terms with definite chiral parity. The chiral even terms are 
proportional to 7^ and 7^75 and are non-zero in the massless quark limit; the chiral odd terms 
are proportional to a^ v and vanish in the massless limit. These cases correspond to the quark and 
antiquark having antiparallel (chiral even) or parallel (chiral odd) helicities. 
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as 1/t 4 and the results do not agree with experimental data 8 . Furthermore, this 
calculation used two-gluon exchange and will give a flat energy behavior. 

Chiral odd components were also included by Enberg, Forshaw, Motyka, and 
Poludniowski together with complete BFKL evolution in Refs. 10 and llj but with 
a different non-perturbative treatment of the chiral odd amplitudes. As we will see 
this gives a ^-dependence of the differential cross section which agrees very well with 
data. Moreover, transversely polarized mesons dominate for moderate momentum 
transfers. (See Section for more details). 

In yet another approach, Ivanov and KirschneJ23 calculated the process using 
a meson wave function with explicit fc-dependence obtained from the photon wave 
function through the use of QCD sum rules. This meson wave function reduces to 
the distribution amplitude when integrated over transverse momentum, and allows 
studying the effects of large quark-antiquark dipoles, since these dipoles have been 
assumed small in other analyses. The calculation used Born level pomeron exchange, 
but additionally included double logarithmic corrections as an approximation of the 
BFKL evolution. Chiral odd contributions were not included. 



2.4. Heavy meson calculations 

Heavy meson production (J/ip and T) is usually considered in the approximation of 
very heavy quarks, in which the quarks have zero transverse momentum and share 
the longitudinal momentum of the meson equally. The first BFKL calculation of 
this process was done by Forshaw and Ryskin.^w ho s olved the BFKL equation ap- 
proximately. Bartels, Forshaw, Lotter and WusthofL then computed the scattering 
amplitude analytically for asymptotically large rapidities using Lipatov's BFKL so- 
lution. This result was compared to HERA data by Forshaw and PoludniowskPSI 
who found good agreement between theory and data for heavy mesons. But these 
authors also considered data on p and <fi production. Surprisingly, the predicted dif- 
ferential cross sections agreed quite well with light meson data, although the heavy 
quark approximation should not be applicable for these mesons. This agreement 
is understood when considering the full BFKL amplitudes'^] as I w iH explain 
below. 

The above calculations were all done in the approximation of asymptotically 
large energies (i.e. rapidities). Technically this amounts to retaining only the leading 
conformal spin in the expansion of the full BFKL amplitude. The conformal spin is 
an integer n, which labels the eigenfunctions and the full amplitude is a sum over all 
n. The only term in the sum that does not vanish in the limit of rapidity y — > oo is 
n = 0, so to simplify the calculations the other terms are often dropped. However, 
it turned ouf " '' I that to describe the cross section for gaps between jets using 
BFKL, it was necessary to keep the subleading terms in the sum. As a consequence, 
Enberg, Motyka, and Poludniowski^ computed the analytical expression for the 

g See Ref.E3for a comparison of the results of Ref.EHto the ZEUS dateP. 
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Fig. 3. BFKL computatiorESl compared to HI data=Ion the differential cross section. The solid 
line shows the LL BFKL result mentioned in the text. This plot is taken from Ref. I^l 

heavy meson production amplitude for all conformal spins. It was found that for 
J 'Ity production the non-zero conformal spins contribute roughly 10% to the cross 
section at the HERA experiments. This calculation has been compared to recent 
data with good agreement, see Fig. EH 

3. General BFKL result for vector meson production 

This section contains a review of the work reported in Refs. I1UI 1111 mentioned 
above. The purpose of these papers was to compute the complete set of amplitudes 
for the vector meson production process jq — > Vq by BFKL pomeron exchange, 
for any momentum transfer i, any meson polarization A, and any vector meson 
V = p^cj),!^, J ftp, T, unifying and generalizing previous work. 

I have already stressed above that due to large logarithms, BFKL evolution of 
the pomeron exchange is needed in order to obtain the correct Regge-type energy 
dependence of the amplitudes 11 . We will also see that the BFKL resummation has 
the effect that amplitudes are always finite, even for zero quark masses — in contrast 
to the Born level two-gluon exchange. Furthermore, the resummed amplitudes are 
dominated by configurations where the pomeron momentum is shared between the 
two gluons, instead of very asymmetrically distributed as for two-gluon exchange. 

In summary, the calculated amplitudes have the features that (i) the pomeron 
exchange has LL BFKL evolution, for all conformal spins; (ii) quark masses are 
included for all vector mesons; (iii) chiral odd photon splittings are considered by 



h See e.g. the comparison with data in Ref. l2l 
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using the constituent quark mass in the perturbative QED wave function; (iv) full 
higher twist distribution amplitudes are taken into account, and both the asymptotic 
and non-asymptotic forms have been tried; and (v) all amplitudes are finite at the 
endpoints (u — ► 0, 1). Thus, the resul ts in this section are very general and contain 
all the earlier calculation d^ l ^ l 'n28 | 29 | 5 ni m various appropriate limits. 

3.1. Impact factors 

We want to compute the helicity amplitudes for the vector meson process, M\\>, 
where A and A' refer to the polarization of the incoming photon and of the outgoing 
vector meson. Since the photon is real, we need only consider transversely polarized 
photons, which we can take to be + polarized. The meson has A' = +, — or 0, which 
means that we need to compute the helicity amplitudes A4+o an d M.-\ 

Because the amplitude is factorized, we may start with the amplitude for two- 
gluon exchange according to Eq. (J2J), and then replace the middle block in Fig. 
by the BFKL pomeron instead of just the two gluons. According to Eq. (J2J), we 
then need the impact factors <?'>'( A )^ V ( A ) and (p q ^ q for the upper and lower blocks 
in Fig. |2 The quark impact factor is a constant, and we choose the normalization 
(p q ^ q — 1, putting all non-trivial factors into the photon-meson impact factor. 

We include both the chiral even and chiral odd photon splittings, writing 
<P = ^ovon + ^odd for all impact factors. The chiral even impact factors were com- 
puted in Refs. IT3l 1241 for massless quarks. Instead of the non-perturbative photon 
wave function of Ref. 1241 we account for the non-perturbative coupling and chiral 
symmetry breaking by calculating the chiral odd impact factors with a perturbative 
QED photon wave function, replacing the current quark mass with its constituent 
mass. Note that all the chiral odd amplitudes vanish for vanishing quark mass. 

Full expressions for the impact factors for all helicities and chiralities are given 
in Ref. 1101 Here, as an example, I will only display one of them, 

4«i Wl ' (+, = K^^eQv /^Hr|*i(H<-|)(e <+> •*<+>•) 

«^/*^ (/ ^ «») + / £*,,(.)) , (4) 

\0 u ) 

where m is the quark mass, and My, fy and Qy are the mass, decay constant 
and effective quark charge of the vector meson, respectively. 4>\\(u) is the twist-2 
distribution amplitude for the vector meson. 

3.2. Solution of the BFKL equation 

The next task is to replace the two-gluon exchange mechanism with the full BFKL 
pomeron. In order to do this, we use the beautiful solution of the BFKL equation by 
Lipatovpl which was obtained drawing on the conformal symmetry of the BFKL 
kernel when expressed in two-dimensional complex gluon coordinate space. The 
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amplitude of a scattering process takes the form of an expansion in the basis of 
conformal eigenfunctions E nM of the BFKL kernel, 

ka V" f a + n 2 /4) exp[(3a s /7r)xn(t/)?/] , ,, u . |)? , 

M = }_ j Jdv [|/a + (n _ 1)2/4][y2 + (n + 1)2/4] (E n ^) ($ 2 \E n , v ) , (5) 

where n is the conformal spin, y = log(s/A) is the rapidity, and Xn{v) is the char- 
acteristic BFKL function given by Xn{v) = 2Re (ip(l) - ij){\/2 + \n\/2 + iv)). The 
factor {$>\\E ntV ) is a schematic way of representing the projection of the impact 
factor on the eigenfunction, which is given by 

1 f d 2 k f 

{$i\E n<v ) = — J ^i(fe) g) J d 2 pid 2 p 2 E n ^(pi,p2)exp(ik-pi+i(q-k)-p2), 

(6) 

where the eigenfunctions are 

/ \ -£+1/2 / / x * x -At+1/2 

^(^2) = (f±=£) ((^) ) , (7) 

with p = n/2 — iu and p, = ~n/2 — iv. Here k and q are transverse two dimensional 
momentum vectors, and p\ and p 2 are space vectors in the complex representation. 

Computing all integrals that can be done analytically, one finds that all six 
hclicity amplitudes take the form 

i ^cvcn/odd ^ \ ^ / 7/7 ^rcvcn/odd/ \ / n \ 

M "' = u 2^ / dv ^ xx' (u) ( 8 ) 



[v 2 + {n- l/2) 2 ][v 2 + (n + 1/2) 2 ] sin(W) 



exp[(3a s /7r)x2n(^)y] r , _ , 
— I a p(v, 2n,q,u;a) 



Here, dd (u) are functions (in which all constants have been absorbed) which 

are combinations of the distribution amplitudes, and which specify the u-dependence 
for different helicity combinations. I a p(v,2n,q,u]a) is given by a line integral, 

I a p(v, n, q,u;a) = - ^{a/2 - C)r(-a/2 - C) t* (i sign (1 - 2u))"^+" 

x (A) [sin7r(a + ^ + C) B(a, p,q* ,u, Q B(/3,Ji,q,u* , C) 



191 

- (-1)" sin7r(a - /* + B K "A*, ?*, «. B(A -ft 9, «*, ] , (9) 

where we have introduced the dimensionless parameter r = 4m 2 /\q\ 2 and the con- 
formal blocks 

r{p + 2 + a + () 



B(a,p,q*,u,0 = (-4«(1 - u ))-(^+2+«+C)/2 2 



r( M + i) 



2i H 2 ' 2 ; ^ + 1; MT^)J- (10) 

The parameters a, /3, and a above take on different constant values for different 
amplitudes (see Ref. 111)1 for a detailed derivation of formulas 1101 ). 
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Fig. 4. Helicity amplitudes differential in u, at fixed y and t. (a) BFKL results compared to the 
massive quark two-gluon result. The two-gluon amplitudes have been multiplied by a factor three, 
(b) BFKL results for three different constituent quark masses, m = My/2 (blue, dashed curves in 
the middle) , and m = My / 4 and m = My . From Ref . 1111 



The nice symmetric structure of these expressions is a consequence of the original 
conformal symmetry of the BFKL kernel. Note that the functions B are not single- 
valued as they stand, but that the complete integrand in Eq. © is. The sum over 
conformal spin includes all terms with even n, although they are rapidly becoming 
very small and only about a dozen are needed for phenomcnological applications. 
However, only the terms with odd n vanish identicall y, in c ontrast to the similar 
case when the vector meson is replaced by a photon!^^ In this case only the 
terms with n = 0, ±2 contribute and the others vanish. The reason for this is not 
known. 



3.3. Properties of the solutions 

The BFKL resummation, in contrast to the two-gluon approximation, makes the 
amplitudes infrared finite in the endpoint region. In Ref. 1^1 the amplitudes 
M^q, and were found to have divergent integrals over u, while in the BFKL 

case these integrals converge. One may think that this regulation of the divergence 
is due to the inclusion of quark masses, but in fact the amplitudes of Ref. 1Q: remain 
convergent in the limit m — > 0; in Appendix B of Ref. 1111 it is shown that the 
amplitudes converge unless one takes the simultaneous limit m — > and y — > Q. 

In Ref.|]2]it is suggested that the endpoint divergences could be suppressed by a 
Sudakov factor limiting radiation off the soft dipole configurations. In Ref.OSIon the 
other hand, it is argued that it is precisely these endpoint divergences that provide 
the solution to the problem that the amplitudes do not scale with t as expected. 
The divergence is taken as a sign that factorization breaks down, and treating this 
correctly provides an extra factor of t. For more details, I refer to this paper. 

The regulation of the endpoint divergences begs the question of how the u- 
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dependence of the amplitudes is affected by BFKL effects. Fig. 0Ji shows the six 
helicity amplitudes differential in u calculated from BFKL and from two-gluon ex- 
change including quark masses, after evaluating the remaining integrals numerically. 
As discussed above the quark masses regulate the endpoint divergence of the two- 
gluon amplitudes, so the divergences are not seen in this plot, but it shows that 
BFKL somewhat suppresses the contributions of asymmetrical dipoles relative to 
the two-gluon curves. 

But how is this related to the choice of constituent quark mass? The physically 
motivated choice for the quark mass is m = My/2, but in our calculation this is 
a non-perturbative parameter which may be adjusted. The answer is that when 
the constituent quark mass is increased the curves are more concentrated towards 
u = 1/2, see Fig.^Jj. This is especially true for the amplitudes that were divergent 
for the massless two-gluon case, .MST 1 , A4°$f, and .M which seem to have some 
memory of said divergences. 

Finally, with the complete BFKL amplitudes we can investigate how the heavy 
quark approximation used earlie J ' l ^v l ^y l w l arises from the more general amplitudes. 
This is studied in some detail in Ref. 11 As was discussed above, the approxima- 
tion consists in taking the quarks to be very heavy and to move together without 
relative motion. This means that their longitudinal momenta are equal, leading to 
the replacement <fi(u) — > S(u — 1/2). Moreover, the transverse momenta of the quark 
and antiquark are zero, which in coordinate space means that the transverse size of 
the pair is zero. These two replacements leave only the .M^^ amplitude non-zero. 
But in fact, it is precisely the amplitude that dominates the cross section 

also for the full BFKL result. The obvious conclusion is that the reason for the phe- 
nomenological succesJ^ of the heavy quark approximation (even for light mesons 
where there is no reason to expect it to work) is due to the selection of the dom- 
inant amplitude in that approximation. This success, however, is limited, because 
there is no helicity flip of the quarks, and the mes ons produced are all transversely 
polarized, at variance with experimental dataP^ 



3.4. Comparison with data 

Now has come the time to consider how the theoretical amplitudes compare to the 
experimental data. In Ref. 1111 we performed a thorough comparison of both the 
BFKL and two-gluon results calculated in Ref. E3 with ZEUS data^ on p and <j> 
production and HI data- on J/ip production. In this section I will therefore only 
highlight a few of the results and conclusions. 

In the leading logarithmic approximation it is not prescribed how to treat the 
running coupling ot s , and the normalization of the cross sections is therefore not 
completely fixed. The running coupling actually appears in two distinct places, in 
the coupling of the two gluons to the impact factors and in the coupling of gluons 
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Fig. 5. Differential cross section for p photoproduction: theory compared to ZEUS dat£0 The 
solid line shows the BFKL result with fixed coupling ctl F = 0.17 and af FKL = 0.25. The dashed 
and dotted lines show the two-gluon result with fixed and running coupling, respectively. 

inside the BFKL ladder. We use a fixed value for the latter, 1 but we use both 
fixed and running a s in the prefactor. The definition of the energy scale in the 
rapidity, y = ln(s/A 2 ), is not prescribed at our level of accuracy either, and we use 
A 2 = My — t. Thus, in the end, we have three parameters in our fits, the two a s 
and the quark mass m, chosen as m = My/2 in comparing to the data. 

The BFKL differential cross sections da/dt are found to be in very good agree- 
ment with the data, see Fig.[5]for p and0for J/ip production. The same parameters 
have been used for all mesons (including </>, which is not shown here). We see that 
J/tp can be reproduced with the full BFKL calculation using the DAs, in addition 
to the earlier calculations in the heavy quark approximation shown in Fig. [3J The 
two-gluon curve with running coupling is in as good agreement as the BFKL curve, 
but with fixed coupling it is clearly too flat. The reason is that the flatter A4+o 
amplitude is much larger in the two-gluon approximation than in the BFKL case. 

This is evident when looking at the measured spin density matrix elements r®j, 
which arc dimensionless combinations of the different M\\> amplitudes. They can 
be measured experimentally through angular distributions of the meson decay prod- 
ucts. This provides the only way of studying the helicity structure of the process; 
to measure individual helicity cross sections one would need polarized beams. 

There are three relevant spin density matrix elements in this process: rgg oc 
|M+o| 2 , oc Re[M* +0 (M ++ - M+-)}, and r 04 n oc Re{M ++ M* + _]. All are nor- 
malized to the sum of all amplitudes. Tqq directly gives the ratio of the Al+o to 
the full amplitude and as seen in Fig. for p production, this is well reproduced by 
BFKL but completely wrong for the two-gluon curve. Thus, even though the cross 

'This is well-motivated and common practice in LL phenomenology. Furthermore, it has been 
arguecPU that the effect of higher order logarithms is to effectively freeze the coupling. 
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Fig. 6. Same as Fig. ISlbut for J/ip production and HI datJ^l 
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Fig. 7. Spin density matrix for p photoproduction: theory compared to ZEUS dateP5 The solid 
line shows the BFKL result and the dashed line shows the two-gluon result. 



section was described by two-gluon exchange, the dynamics is not correct and the 
dominance of transverse mesons seen in data is not reproduced. 

However, we also find that r_ xl is somewhat too negative for the BFKL curve, 
while the real failure shows up for r^: both the BFKL and two-gluon curves have 
the wrong sign! This shows that the M+o amplitude has the wrong sign in both 
calculations, and consequently that it is not under control. There are large end 
point contributions in the calculation, which may, as in Ref. signal that the 
used factorization is breaking down. 

An interesting observation, however, is that if the quark mass is increased to 
double its physically motivated value, then the sign of tJq is switched to the cor- 
rect one (the cross section is much less sensitive). This is illustrated in Fig- El This 
increase of the mass may be somewhat artificial, but it serves to show an interest- 
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ing point: Increasing the quark mass effectively leads to a stronger suppression of 
large dipoles (through the Bessel functions in the impact factors). There are several 
possible sources of such a suppression: saturation effects, Sudakov suppression of 
radiation off the quark and antiquark, or inclusion of the size of the vector meson 
in the wave function such as in Ref. 1271 

For J/tj} production all r^ 4 are consistent with zero, but the errors are quite 
large. Within the heavy quark approximation all r° 4 identically vanish, since only 
the amplitude is non-zero. However, the full BFKL calculations predict non- 

vanishing rf 4 , and these results are also compatible with the J/tp data. 



4. Conclusions and outlook 

In summary, the BFKL calculations reproduce the differential cross sections for all 
vector mesons very well, and also the dominance of transversely polarized mesons 
and relative importance of the different helicity amplitudes. But they do not cor- 
rectly account for the sign of the amplitude for longitudinal mesons. The two-gluon 
calculations, on the other hand, predict a much too large fraction of longitudinally 
polarized mesons but can despite this describe the differential cross section. How- 
ever, the problem of endpoint divergences remains. 

More theoretical work is clearly needed to understand the process, but in any 
case it has been demonstrated that BFKL is an important ingredient in the calcula- 
tions and, importantly, makes the amplitudes finite. The results have also confirmed 
that the chiral odd components of the photon wave function are important to de- 
scribe the data. It seems, however, that we may require a larger suppression of 
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large dipoles, and the role of the factorization into a hard process and soft vector 
meson formation needs to be investigated. Computing the corresponding processes 
in electroproduction would also give more constraints on the helicity structure. 

On the experimental side, it would be helpful to have data for larger momentum 
transfer and with better statistics. If, with more precise data on J/ip production, 
one would see deviations from zero in the spin density matrix elements, it would be 
a clear sign that the heavy quark approximation is not appropriate and that one 
should instead use the full framework discussed here. 
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